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Abstract 

We consider a porous solid covered with a water film (or with a drop) 
in situations where the liquid is pumped in, either spontaneously (if the 
porous medium is hydrophilic) or mechanically (by an external pump). 
The dynamics of dewetting is then strongly modified. We analyse a few 
major examples a) horizontal films, which break at a certain critical 
thickness b) the "modified Landau-Levich problem" where a porous 
plate moves up from a bath and carries a film: aspiration towards the 
plate limits the height H reached by the film c) certain situation where 
the hysteresis of contact angles is important. 

Keywords: porous media, wetting, capillarity, contact lines. 

PACS: 68.45. Gd, 68.15.+e, 47.f5.-x. 

1 General aims 

The current picture of dewetting processes on flat, smooth surfaces is relatively 
clear 0, although inertial regimes still raise some questions 0. But if the 
surface is porous, the situation is very different. A number of practical problems 
belong to this class: 

(i) in the 4-color offset printing, a sheet of paper may be at one moment 
covered by a water film. The water must be removed very fast when the sheet 
reaches the desired printing roller: this dewetting process is assisted by a spon- 
taneous suction of the water into the hydrophilic paper || . 

(ii) many industrial processes involve the rapid passage of a solid film, into a 
liquid bath. When the film leaves the bath, it drags a Landau-Levich film |4|, || . 
However, if the film is porous, and sucks in the fluid, there is a height H at 
which the Landau-Levich film disappears: what is the value of HI Ejj. 

(iii) when we paint with a brush, we are in fact transferring a liquid from a 
porous medium onto a flat surface - the reverse process. 
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Our aim, in the present text is not to cover all these complex systems, 
but to choose a few model examples where the effects of aspiration (towards 
the solid) are non trivial. The essential parameter describing the aspiration is 
the normal current J (volume/unit area/unit time). When we pump through a 
given thickness of porous solid, with a prescribed pressure drop, J is controlled 
by Darcy's law. When we have a spontaneous suction, over a time interval t, J 
is related to the Washburn equation the scaling structure of J (for a pore 
volume fraction ^1/2 and wetting angles ~ 1 radian) is simply 



J^j) . (1) 

In eq. (Qj), the parameter D has the dimensions of a diffusion coefficient: 



D = -d=V*d, (2) 
V 

where d is the pore diameter, 7 is the surface tension, r\ the viscosity of the 
fluid, and V* — j/rj is a characteristic capillary velocity. In all the model 
systems to be discussed below, we shall assume that J is prescribed and is time 
independent. (One case where the time dependence may matter is discussed in 
rcf. j(|). We also assume that the pore diameter d is very small: then we may 
describe the capillary hydrodynamics in a continuous picture. 

One natural feature of a porous surface is to pin down a contact line - i.e. 
to bring in a certain hysteresis in the contact angles. If the equilibrium value of 
this angle (as given by Young's law) is 9 e , the receding angle is 6 r < 9 e . Two 
distinct cases can be found in practice: 

a) 'strong' pinning 6 r = 0. 

b) 'weak' pinning 9 r > 0. 

In this paper, we start (sections |^ and ^ by a discussion of pumping effects 
on surfaces with no hysteresis^: this is not usual, but conceptually important. 
We then proceed to some cases of hysteresis inducing weak or strong pinning. 

In section ^J, we consider horizontal films, which thin out and become unsta- 
ble at a certain critical thickness. Then we proceed in section ^| towards vertical 
plates pulled out of a liquid: the modified Landau-Levich problem. In section [| 
we allow for hysteresis, first solving the macroscopic problem of a pinned drop, 
then reconsidering the vertical plate situation. Section^ discusses some possible 
extensions and some limitations of the present work. 



2 Horizontal films 



1 In order to have a good determination of the thermodynamic angle 9 e , we do not neces- 
sarily need an ideal surface without defects. As shown by J. F. Joanny and P. G. de Gennes 
[J. Chem. Phys., 81, 552 (1984)] in the case of "regular" defects, we need only a surface with 
irregularities below a certain strength threshold. In sections M and pi the local dissipation 
which might be associated to these defects is assumed to be negligible. 
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2.1 A macroscopic puddle under aspiration 

This case is described on fig. [|. It has been studied recently theoretically and 
experimentally by Bacri and Brochard ||. Here we present only the case with 
no hysteresis, because it provides a good introduction. 

Let us denote by S the spreading parameter jso — (7 + Isl) where 7 is 
the liquid surface tension, and jsl (resp. jso) the interfacial tension of the wet 
(resp. dry) solid. We assume S to be negative (partial wetting). The equilibrium 
puddle has a thickness eo which results from a balance between gravity (tending 
to decrease eo) and capillarity (struggling for a minimal exposed area). For a 
large puddle (radius R ^> eo) on a partially wettable substrate, this corresponds 
to the following energy: 



, 1 2 O , 

/ = ^pge — + — (7 + Isl - lso) (3) 

where fi is the liquid volume, Q/e the contact area, p the density, and g the 
gravitational acceleration. Minimizing (||) with respect to e gives a classical 
formula for eo (which is typically of order 1mm) : 




e = W2— . (4) 



Let us now pump on this structure, imposing a small current J (and as- 
suming zero hysteresis). Then e remains close to eo. The horizontal shrinkage 
velocity U is related to J by the conservation equation: 



2TrRe Q U = nR 2 J, (5) 

from which we deduce 

All the structure is slightly distorted. For instance the dynamic contact angle 
9d is slightly smaller than 9 e , because there is a Poiseuille flow in the vicinity 
of the contact line. The necessity for this flow can be perceived starting from 
the opposite assumption: if the liquid in the puddle was flowing downwards, 
uniformly with velocity J, the lateral velocity, deduced from the position of the 
contact line, would be U — J/9 e ^ U . We can estimate the difference 9 e — 9^ 
from a standard dissipation argument jlj] : 

TE = ^-(U - U) 2 log — = (C/-C/)7(cos0 d -cos0 e ) (7) 
9d a 

(where we have assumed 9 e and 9d to be small, since this is the most universal 
case). The argument in the logarithm is the size ~ eo of the Poiseuille region 
(near the rim), divided by a molecular length. 
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2.2 Nanoscopic pancakes or films 

We now pump on a very thin film (e = nanometers), and assume first that the 
solid surface is partly wettable (i.e. jso < 7sl + l)- At these small scales, the 
film energy F(e) per unit area can be written as fllPf 

F(e)= 1S L+l + P(e), (8) 

where, for e much larger than the molecular size ao, -P(e) is controlled by long- 
range Van der Waals forcesf]: 

P ( e ) = TT-^ (e»a ). (9) 

We will assume that the Hamaker constant A is positive. For shorter distances 
(e ^ ao), the precise form of P(e) is generally unknown, but in the limit e — > 0, 
we have P(e — > 0) = 750 — 7 — 7s l since we must recover F = jso (dry solid). 
The general aspect of F(e) for partial wetting is shown on fig. ^a. 

Our film is locally stable when the curvature F"(e) is positive jy]: thus, 
in the case of fig. ||a, we can get down without changes to a certain thickness 
e s , such that F"(e s ) = 0. But, when e reaches e s , a spinodal decomposition 
is expected. If we stopped the pumping at this moment, this would ultimately 
lead to a coexistence between dry regions and finite, macroscopic droplets. But 
if we persist in pumping, the early droplets (of thickness ~ e s ) will fade out 
before any major coalescence process. 

A similar discussion can be given for the case of complete wetting ("/so > 
Isl + 7). The corresponding plot of F(e) is described in fig. ||b. In the absence 
of pumping, the liquid forms a "pancake" of nanoscopic thickness e = e p . The 
equilibrium condition determining e p is || 

e p n(e p ) + P(e p ) = S (10) 

where life) = -dP/de. Eq. @ may be interpreted as a balance of forces on 
the contact line: the Young force S = jso — (jsl + 7) is equilibrated by the 
pressure due to Van der Waals attractions. (The graphical construction of e p is 
shown on fig. ||b: e p is the point for which the tangent to F(e) intercepts the 
vertical axis at 750-) 

As a consequence of this, when pumped, a continuous film would tend, at 
e < e p , to break into a mixture of such pancakes and dry regions rather than 
thin out uniformly. Each pancake will then shrink by a process reminiscent of 
the puddle [e.g. eq. (||)]. 

Of course, all this discussion of nanoscopic objects is rather unrealistic, since 
most porous media have pores larger than the structure described here: our 
continuum description can be adequate only if our porous surface was based 
on an extreme permeation membrane, with pore diameters of a few angstr0ms! 
But, in spite of this strong limitation, the discussion is (we think) conceptually 
helpful. 

2 We only retain non-retarded Van der Waals interactions: the retarded regime occurs for 
distances larger than those that will be considered in this paper. 
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3 The modified Landau-Levich problem 



3.1 A reminder on non porous systems 

We pull out a vertical plate from a bath (fig. ||a) at velocity V. If the plate 
is partially wettable (equilibrium contact angle 6 e ), and the velocity V is very 
low, we retain a dynamic contact angle 9d close to 9 e . If we increase the velocity 
V, we have a non trivial relation V(9d). Of course V{9d) vanishes at 9d = 9 e , 
but V also vanishes at 9d — > 0, because the dissipation in a thin wedge is very 
large. The simplest (crudest) form, valid for small angles, is Q 

V{9 d ) = const. V* {f e - 9 2 d ) 9 d . (11) 

This has a maximum: 

V = V m = const. V*9 3 e . (12) 

If we impose a velocity larger than V m (as usual in fast technical processes) 
we cannot retain a contact line: the plate drags with it a film of finite thickness 
ex,. This thickness was computed in a classic paper by Landau and Levich Q 
(for viscous regimes): 

V \ 2 /3 



eL = b {v^) (13) 

where b is related to the curvature of the underlying meniscus (6 ~ 1 mm). The 
height of the adjustment region (beyond which the film becomes uniform) is of 
order biV/V*) 1 ' 3 . 

The above results for the film thickness and the adjustment region are also 
valid for a completely wettable plate. In that case, however, there is no velocity 
threshold and the plate drags a film as soon as V > 0. 

Eq. has been amply verified, and extended up to inertial regimes ||. 
Our aim here is to see how the film is modified under pumping, in a viscous 
regime of small capillary numbers (Ca = V/V* <C 1). 



3.2 Pull out of a surface under pumping: Macroscopic 
flow 

We now consider the pull-out of a substrate able to suck liquid in as it drags 
a film from the bath. This modified Landau-Levich problem was discussed in a 
simplified form in ref. Q, and, for completeness, we start by a brief summary 
of the main results. 

Let us assume that we impose a uniform pumping current J on the film, 
J being small compared to the drag velocity V. The aspect of the flow is 
pictured on fig. [|b: as in the usual Landau-Levich film, there is a meniscus 
from the bath and an adjustment region leading to a thickness e/,. But then 
on, aspiration comes into play, and we expect that most of the film will have a 
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simple hydro-dynamic flow (plug flow) in the external film, without any pressure 
gradients, where the flow lines are just tilted by the angle 6 = J/V <Cl. As a 
consequence, the Landau-Levich film thins out and ultimately disappears at a 
finite height H — el/0. 

However the validity of this non-dissipative solution is necessarily restricted 
to macroscopic scales, and is not acceptable very close to the contact line, where 
one must take into account local wetting properties such as the spreading pa- 
rameter S: the slope of the profile must change. It is our aim in the rest of 
this section to understand the structure of the film in the vicinity of the contact 
line, both in the case of complete (S > 0) and partial wetting (S < 0). 



3.3 Hydrodynamic equation and boundary conditions 

Let us denote z the vertical distance measured downwards from the contact line. 
In the lubrication approximation, we can describe the velocity in the film (in 
the z-direction) by the following parabolic form: 

v z {x) = -V + 4(z 2 - 2ex), (14) 

dictated by the no-slip condition on the plate, and the absence of tangential 
stress at the free surface (note that in this formula, V\ and e are functions of 

*)■ 

For a steady state situation, conservation imposes that the upwards flow 
rate — J Q v z (x)dx at height z balances the amount Jz of liquid that disappears 
above due to the pumping: 

J z = -J o v z (x)dx =e(v+ ^V^j . (15) 

The Poiseuille flow ( |l4| ) is driven by a pressure gradient dp/dz, where p contains 
contributions from capillarity and Van der Waals forces 

V = -^-^ (16) 

and A is a Hamaker constant. 
Writing 

dp = ? d 2 v z = 2??yi 
dz dx 2 e 2 ' 

equation (|l^) can be rewritten as 



9z _ e 2 [ d 3 e d_ f ( 
where, as before, 8 = J/V is a small angle (# <C 1 
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Equation ( JL8[ ) is reminiscent of the equation discussed in ref. m for precursor 
films over a wettable surface, and in ref. p"2| for an advancing liquid on a partly 
wettable surface. But there are two differences: a) pumping introduces a new 
term (—8z/e) on the left hand side, b) the direction of the driving velocity V is 
reversed (in the frame of the plate, the liquid front moves back). 

In order to single out the physical solution of eq. (|l8|), it is necessary to 
impose the three following boundary conditions: 

i) e(z) = at z = ; (19) 

ii) lim — — ; .„ - ie'(z) = - ; (20) 
' z^o 127re(z) 2 2 w 7 ' y ' 

Hi) e reaches a linear profile when z goes to infinity. (21) 

The boundary condition (ii) reflects the local balance of forces on the contact 
line ||. Since it involves the spreading parameter S, it will introduce a dis- 
tinction between situations of total wetting (S > 0) and situations of partial 
wetting (S < 0). 

Far from the contact line, as expected, the profile merges into the macro- 
scopic profile of angle 9 described in the previous section: in eq. dl8|), for large 
thicknesses and vanishing curvature, e ~ 9z + C is easily seen to be the asymp- 
totic profile (C a constant). But a more thorough analysis of the asymptotic 
behaviour of eq. ( |l8| ) |l3| shows that any profile of angle 9 is not convenient, and 
that we must necessarily set C — 0. Hence, remarkably enough, if the asymp- 
totic profile were known with sufficient accuracy, by a straight extrapolation of 
it one could guess the actual microscopic position of the contact line, regardless 
of any complicated local behaviour. 

Bearing this essential fact in mind, we now turn to solving eq. ( |l8| ) in the 
case of a completely wettable substrate. 

3.4 Wettable surface 

We first consider the modified Landau-Levich problem on a wettable surface. 
From now on, to make later explanations clearer, the two terms of the l.h.s. 
of eq. ( |l8| ) will be respectively referred to as the 'dynamic term' (since the 
presence of the "1" is directly due to the pull-out) and the 'aspiration term', 
and the terms on the r.h.s. respectively as 'capillary' and 'VW' (Van der Waals) 
terms. 

We start from the 'far-field' region and look for characteristic lengths to 



scale eq. (18). In these scales, we need both the dynamic and aspiration terms 
to be of order one (to give rise to the asymptotic profile of angle 9) , plus a third 
term to describe the departure from it. One can show that the only consistent 
choice is to keep the VW contribution of order one. This imposes the following 
scalings: 

Z\ = z/ni with fix = £ (22) 
ej = e/Ai with Ai = 91 (23) 
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where £ is a crucial length of the problem defined by 



pr* \ 1/2 

T 



(24) 



and a an atomic length (of a few A) defined as a 2 = ^A-. We note that I, 
though much greater than a, remains in the microscopic domain (J/V* = 1 0~ 4 
gives a few hundred A). 

Our central equation (|l~8|) rewrites in these scaled variables 



ei 



(25) 



where 



1 ^*<? 3 

£ =3^ (26) 

is a small parameter of the problem^: e«l. We thus see that in these scales, 
capillarity appears as a perturbative term to a behaviour ruled mainly by the 
pull-out velocity aspiration and VW. However, we immediately notice that the 
capillary term, though of order e, contains the highest derivative in the equation, 
implying that we are in presence of a singular perturbation problem. 

Such problems have the interesting property that even for arbitrary small e, 
there always is a region, called inner region or boundary layer, where the true 
solution dramatically differs from the simple e = solution^ (usually named 
outer solution). For the interested reader, extensive coverage of this subject can 
be found in refs. [l5| [l6|, |l7). In the following presentation, we will try to 
emphasize physical content rather than enter into technical developpements ||l3f . 

In our problem, the inner region stands near the contact line, at z = 0, and 
capillarity is called to play an important role there. The outer region, on the 
other hand, extends from the limit of the inner region (crudely) to infinity. 

The (first-order) outer solution is looked for by setting e = in eq. (p5|). We 
obtain the exact form 

6l '° ut= K + Q evy^dy 1 (27) 

where K is a constant that will be determined later by matching with the in- 
ner solution. One can check that this outer solution respects the asymptotic 
boundary condition (|2l|). 

In the inner region, the essence of singular perturbation techniques is to 
rescale the independent variable with the help of a power of e in order to change 
the dominant balance of the equation. Guided by the form of the boundary 

3 To put numbers, taking = 10 — 2 , V = 0.7 m.s - 1 and V* = 70 m.s - 1 , gives e ~ 3.10 - 5 . 

4 Or even from any perturbative expansion in the powers of e. 
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condition ( ppf ) which balances capillarity and VW at the contact line, we choose 
the following new variables: 

Z2 = z/fj, 2 with fi 2 = e 1/2 £ (28) 
e 2 =e x A 2 = Ai . (29) 

The expression of ^ 2 reveals that the boundary-layer near the contact line ex- 
tends over a distance of order e 1 / 2 £. 

In these inner variables, eq. Hlq) takes the inner form 



e 1/2 - e ?=^' 2 "-§, (30) 
e 2 e 2 

where the dynamic term and the aspiration term remain small, and the be- 
haviour is dictated mainly by a balance between capillarity and VW. Solving to 
first order (by setting e = 0, and integrating twice), we get z 2 as a function of 
e 2 : 

Ce 3 + 2D e 2 + \ 

The lower limit of integration in ( |3l| ) is set to by the boundary condition (|l9|) , 
while application of (20) brings 

(negative for our wettable surface). 

We now proceed to the asymptotic matching jl4|, |l6|, [Tt]] of the inner and 
the outer solution. The first order matching condition writes 

lim e 2 (z 2 ) = lim e 1 (z 1 ). (33) 
The l.h.s. admits a finite limit only if C takes the specific value 

C= 2 -^{-Df'\ (34) 

The value of this limit is then e p /8£, where e p is the thickness of the static 
pancake [eq. (|lo|)l that the liquid would form at rest on the surface, without 
aspiration. We then equate this value with the r.h.s. limit of the matching 
condition, which is equal to 1/K, and find: 

K= 6 1. (35) 

The solutions found for the inner [eq. (511)] and outer [cq. (p7[)] regions, 
together with eqs. (32), (34) and (pq), give us a complete description of the 
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profile of the film over a wettable substrate. Figures |]a and [|b present the 
profile that is obtained, respectively in the close vicinity of the contact line, and 
farther away in the outer region. 

Our results call for a few important comments: 

i) When plotted in orthonormal coordinates, the liquid takes the shape of 
a long, almost flat, strip at the approach of the contact line. This is naturally 
a remembrance of the usual precursor film that is encountered in the wetting 
process of a (wettable) surface by a drop ||. In physical units, the thickness of 
this tongue is indeed similar to that of a precursor, both being of order e p ^\ 

ii) We are now able to give an order of magnitude for the distance over which 



the flow merges with the simple, non-dissipative, flow of sec. 3.2. This merging 



is described by the outer solution e i, out , and in fig. ^b we see that it occurs over 
a distance of a few ^'s [eq. , that is to say on microscopic scales. Practically, 
only the macroscopic tilted flow of angle 9 is observable, thus validating the 
simple analysis of ref. 

iii) One question remains: what determines the actual position of the contact 
line ? We know from the previous section that a delicate hydrodynamic balance 
entails that the contact line necessarily lies on the position obtained by a straight 
extrapolation of the asymptotic linear profile of angle 9. If this asymptotic line is 
shifted, then the whole profile is shifted by the same amount, without change: all 
the positions of the contact line are equivalent. But this degeneracy is naturally 
broken when we attempt to match the asymptotic profile with the meniscus 
originating in the bath from which the plate is pulled out. This matching, by 
fixing the asymptotic line, also determines the contact line position. 



3.5 Partially wettable surface 

We now consider the situation of a partially wettable surface. In fact, physical 
intuition encourages us to think that this will not be a mere extension of the 
wetting case. In the pull-out process, due to aspiration, the contact line is 
forced to recede on the plate. For a wettable substrate, we have shown that this 
dewetting occurs with the liquid leaving, just after the contact line, a thin strip 
similar to a precursor film. On a partly wettable surface, however, the liquid is 
not inclined to leave any film at all, but rather to adopt a sharp profile which 
is in obvious contradiction with the asymptotic line of small angle 9 <C 1. How 
does the system resolve the dilemma ? 

We analyse this situation along the same lines as previously. In the (zi,ei) 
set of variables [eqs. (|2^) and capillarity is still a singular perturbative 

term, but is called to be important near the contact line. The solution will con- 
sequently still exhibit a boundary-layer structure: the outer solution [eq. (p7|)] 

5 The appearance of a static quantity as e p is not a surprise: in the inner region, the 
dynamic term was totally ignored. If we increase the dynamic drive (by increasing V), this 
static approximation can still be made, but the region over which it remains valid (~ e 1/l2 £) 
diminishes accordingly. 

6 In the limit J — > 0, the merging length (~ i) increases towards infinity, but one must 
notice that at the same time 9 — > at a faster rate: there is no contact line anymore, and we 
go back to a classical infinite Landau-Levich film. 
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must be abandoned near in the vicinity of the contact line, where the inner 
solution given by eq. (|3l|) takes place. 



However, in eq. (|3lD, an important change occurs: the constant D = — \¥y?-~ 
[eq. (|3^)] is now positive, because the spreading parameter S is negative in the 
partial wetting case. We then go on and try to match the inner and the outer 
solution with the condition (|33|), but this appears to be impossible: whatever 
the value of the constant C in the inner solution, the l.h.s. of eq. (|3^) does not 
admit any finite limit ^ 

This mismatch takes its roots in the discrepancy between the behaviour of 
the profile in the inner region, where it is similar to the steep, hyperbolic profile 
of a static drop on a partially wettable surface |)| , and the outer region where 
the profile tends to join the low-angle asymptote. Therefore, a third region, 
called intermediate region, must appear between the inner region and the outer 
region, so as to bridge the gap. 

For the intermediate region to play its role, one is led to keep capillarity, 
VW and the dynamic term of order one. To achieve this, both the independent 
and the dependent variables must be scaled with appropriate powers of e. We 
find: 

z z = zj V L Z with /i 3 = 3" 1/6 e 1/6 /ii (36) 
e 3 = e/A 3 with A 3 =3 1/6 e- 1 / 6 Ai, (37) 

so that the 'intermediate equation' takes the form 



e 3 



Deriving an analytical closed-form solution of this strongly non-linear equa- 
tion reveals to be a difficult task. However, a satisfying approximate solution 
can be constructed from a linearized form of eq. (|3^), and be matched through 
a numerical procedure with both the outer and inner solutions to determine all 
the constants that remained undetermined (yielding values of C and K that are 
close to 0). f\ As can be seen on fig. ||, the profile now presents a 'bump' that 
allows a transition between the two extremal regions. 

Let us focus our attention on the evolution of this bump when e becomes 
smaller, for instance when increasing the pull-out velocity V. The length 
and thickness of the bump are respectively of order and 
A 3 = 3 1 / 6 e _1 / 6 A 1 . Therefore, compared to the characteristic lengths /ii and 
Ai in the outer region, the bump's length shrinks, while, as is shown in fig. ^|, 
it gets steeper because of the negative power of e appearing in A 3 . Such a 
steepening appears coherent as an attempt to reconcile conflicting evolutions: 
when V is incremented, the asymptotic profile flattens more and more (9 — *■ 0), 

7 It can he shown that even with the help of more refined matching procedures than that 
given by eq. (B3h, no matching can be performed between the inner and the outer solution. 

8 For the sake of conciseness, a technical presentation of the procedure ]l3| , unessential to 
understand the physical picture, has been omitted here. 
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in increasing contradiction with the 'natural' tendency of the profile to display 
strong slopes in the inner region. 

We conclude by pointing out that, as in the wetting case, the merging dis- 
tance with the macroscopic flow of angle 9 is of order I: for practical purposes, 
there is no visible deviation from the uniformly tilted flow. Indeed, it is an 
interesting fact that, except for a small region near the contact line where the 
profile is submitted to the influence of wetting properties (as the spreading pa- 
rameter), most of the flow is exactly the same in the wetting and the partial 
wetting case. This is of course to be imputed to the far off-equilibrium issue 
that we are dealing with: aspiration by the plate forces a contact line to appear 
in situations where the pull-out is so strong that there should not exist any. 

4 Effects of hysteresis 

4.1 Pinned, macroscopic, horizontal drop 

Let us consider a macroscopic drop on top of the horizontal upper surface of a 
porous medium with hysteresis. The drop is submitted to a uniform downwards 
aspiration (current J) while its contact line is strongly pinned (9 r = 0) at x — 
(the x-axis is along the horizontal surface). Initially, the (macroscopic) profile 
near the contact line, e, is wedge-like (with a contact angle 9o)- e(x,t — 0) = 
9oX. How does the profile evolve under the combined effect of aspiration (which 
would force the contact line to recede) and pinning (which prevents it from 
doing so) ? 

In the lubrication approximation, the horizontal velocity field inside the drop 
writes, using a vertical z-axis: 



This Poiseuille flow is driven by a pressure gradient dp/dx where p = 
—^d 2 e/dx 2 (for the macroscopic drop under consideration, Van der Waals forces 
can be ignored, except in a very small region near the contact line) . 

Writing dp/dx = r\d 2 vjdz 2 = 2rjU/e 2 , we arrive at 




(39) 




d 3 e „ U 



(40) 



The horizontal flow rate, Q = Hf v(z) dz, is therefore given by 



Q 




(41) 



Inserting eq. ( j4l| ) into the conservation equation 




(42) 
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we find the evolution equation of the profile: 

\- J -I —V p 4 = 



at ' J + ^U FV ^J=°- ™ 



For a small pumping current J, e remains close to 9qx and eq. ([43j) can be 
linearized: 



(one can check that the above simplification is actually valid for J <C V™#o)- 

Far away from the contact line, we expect that the profile is not significantly 
different from that obtained in the absence of pinning, which would merely be 
a global downwards translation of the initial wedge by a distance Jt (yielding 
e = 6qx — Jt)- It is therefore natural to look for a solution to eq. (f44| ) of the 
form 

e(x,i) = 9 x-JtF(x,t), (45) 

where the function F describes the effect of the pinning, and is such that, at any 
given time t, F(x,t) ~ 1 for x large enough. We look for a similarity solution 
F(x,t) = f(u), with a variable u of the form x/t a . Replacement in eq. ( fH| ) 
imposes that 

u = 3—^5- , (46) 
V*6lt v ' 

and that / verifies 

-/(«) + «/'(«) + 1 = ^ ( U 3 /'») . (47) 

[The constants in eq. ( |46| ) have been chosen so as to cancel out all physical 
quantities from eq. (§^j)]. 

Since we require that for t > 0, e = at x — 0, we have the boundary 
condition f(u = 0) = 0. Also, for the reasons exposed above, far away from the 
contact line, we must have f(u — * +oo) = 1. We finally impose zero horizontal 
flow at x — and x — * +oo. 

The behavior of / is shown on fig. [Tp,. For u — > 0, / is not analytical. For u — > 
+oo, / presents decreasing oscillations, and behaves like a linear combination of 
e~i ul/3 cos and e~t" 1/3 sin . The essential feature of / is 

that it approximately reaches its asymptotic value (/ = 1) for u of order unity. 

Therefore, we can conclude that the effect of the pinning is felt only in a 
region of breadth £(i) = ^V*6^t. and that beyond this "healing length", the 
drop recovers the profile it would have adopted in the absence of contact line 
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pinning. Finally, the initial (macroscopic) contact angle of the drop is slightly 
reduced during the pumping: 

M J 

T^Wes < 48 > 

[justifying in retrospect the linearization of eq. flls))]. The global shape of the 
drop is shown on fig. ^b. 

4.2 Pull out with hysteresis 

We now return to the modified Landau-Levich problem with the geometry of 
fig. ||b and we allow for some hysteresis in the contact angle. 

We start with the case of strong pinning, where the receding angle 9 r is 
equal to zero. Let us denote F(V) the pinning force felt by the contact line 
at velocity V. Strong pinning implies that, for a flat macroscopic film, the 
force required to depin the contact line, equal to F(V — 0), is larger than the 
capillary force -fsL + 7- P Formally, this is exactly similar to what we had in 
our discussion of section B , with a wettable surface displaying no hysteresis and 
ISO > 7SL + 7- The spreading parameter S should simply be 'renormalized' to 
S = F(V) — ~fsL + 7- hi the region close to the contact line, the liquid forms a 
thin strip, with a thickness of the order of the pancake thickness e p . Here, e p is 
defined by the force balance 

F(V)-(j SL + 1 ) = P(e p ) + e p U(e p ). (49) 

We do not know much about the velocity dependence of the pinning force F(V) 
at large V. For simplicity, we shall assume that F(V) is an increasing function 
of V. Thus F(V) > F(0), and the left hand side of eq. ( |49| ) is constantly positive. 
Hence, eq. ( ff9[ ) admits, in principle, a solution e p (V). 

Again e p is very small, and we expect the visible profile to be dominated 
by the macroscopic solution: dissipation occurs only within a finite vertical 
distance I. At greater distances, we return to a uniform flow of tilt angle 6. 

Let us now turn to the case of weak pinning (6 r > 0). Here, our pinning 
force F(V) is weaker at low V: 

F(0) - (isl + 7) < 0. (50) 

But we are dealing with relatively high velocities V at which F(V) may become 
larger. We are thus led to two possibilities: 

a) if F(V) > 7 + 7sl, we return to the case of strong pinning. 

b) if F(V) < 7 + 7sl, we expect a profile similar to what we discussed 
in section || for partial wetting: at the approach of the contact line, the profile 
thickens, forming a bump. 

In both cases, again, the transition to the macroscopic flow occurs on a 
microscopic region of size I. 

9 In our notation, the capillary force "fso na s been absorbed in F(V). 
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5 Concluding remarks 



We conclude with the following remarks: 

1) A number of interesting experiments are conceivable with a pumped 
porous medium: for instance, in pull out, measuring the finite height H achieved 
by a Landau-Levich film. The surprise is that H is dictated by the uniform tilted 
flow: 

B*%*.^, (51) 

and is insensitive to the delicate dissipation processes taking place near the 
triple line: the "near" field is of nanoscopic size for most cases. And in practice, 
the finite size of the pores may provide a more important cut off than the length 
i discussed in sections (Q) and (||). 

2) There are interesting questions related to the pinning forces F(V) at large 
velocities. Here, the contact line is not necessarily the end point of a liquid 
wedge: it may represent the border between a dry region and a thin pancake. 
But the application of these ideas to porous media, with pores of finite diameter, 
seems remote. 

3) If we want to avoid these difficulties, we can think of pumping a liquid 
droplet, on a non porous surface, by evaporation: this was indeed the starting 
point of a series of experiments by the Chicago group |Q . We benefit by having 
a smooth solid surface. But evaporation has its own complications: temperature 
gradients induce Marangoni flows, which are complex and not very instructive. 

4) Returning to the pull out problem, we should also insist on another limi- 
tation of our discussion: we restricted ourselves to steady state regimes. At high 
velocities, the liquid wedge might decide to emit a periodic train of droplets (or 
pancakes): these droplets would then be pumped out by the current J. A tap 
delivering a low output of water is a good example of periodic droplet emission. 
This idea is attractive. However, because of the nanoscopic size of the objects 
involved, the effects may be hard to observe. 
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Figure 1: Macroscopic puddle under aspiration. 



Figure 2: Film energy versus film thickness in the case of: a) partial wetting, 
b) complete wetting. 



Figure 3: Pull out at velocity V of a vertical plate from a bath: a) The classical 
Landau-Lcvich problem, b) the modified problem, with a porous plate inducing 
an aspiration current J. Most of the film then displays a tilted profile of angle 
9 = J/V, with a plug flow velocity field, and terminates at a height H = e L /9. 
In the present article, we investigate the microscopic structure of the film near 
the contact line (dotted circle). 



Figure 4: Film profile in the complete wetting case, for e = 2.3 1CT 5 and 
S/j — 0.2 (length unit is t on both axes), (a) Vicinity of the contact line, (b) 
Outer region. After a few £'s, the film merges with the macroscopic tilted flow 
of angle 9 (asymptote). Note, in each plot, the very different scales between 
the two axis : the one giving the distance from the contact line should be 
considerably stretched to get a representative idea of the actual film profile. 



Figure 5: Film profile near the contact line in the partial wetting case, for 
9 = 10~ 3 and e = 2.3 10~ 5 (length unit is I on both axes). The profile presents a 
'bump' corresponding to the intermediate region. Farther away from the contact 
line, the profile approaches the asymptote in a way similar to the wetting case 
(not shown). 



Figure 6: Bump steepening in the film profile (partial wetting) as the pull-out 
velocity is increased (decreasing values of e). 



Figure 7: Macroscopic, pinned, drop under aspiration, (a) Shape of the function 
f(u) (see text), (b) Deformation of the liquid wedge near the line. The grey 
line represents the initial profile. The effect of pinning is felt inside a region 
of width £(i); beyond, the profile (solid line) "heals" and recovers the shape it 
would have had in the absence of pinning (dotted line). 
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